Active vs passive scalar turbulence 
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Active and passive scalars transported by an incompressible two-dimensional conductive fluid are 
investigated. It is shown that a passive scalar displays a direct cascade towards the small scales 
while the active magnetic potential builds up large-scale structures in an inverse cascade process. 
Correlations between scalar input and particle trajectories are found to be responsible for those 
dramatic differences as well as for the behavior of dissipative anomalies. 
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The transport of scalar fields by turbulent flows is a 
common physical phenomenon. The dynamics of the ad- 
vecting velocity field often depends on the transported 
field, which is then dubbed active. That is, for instance, 
the case of temperature, affecting the velocity via buoy- 
ancy forces. Conversely, situations where the flow is not 
influenced by the scalar, e.g. for the concentration of 
a dilute tracer, are referred to as passive. Progress has 
been recently made for the passive case (see Ref. [1] and 
references therein). These studies put renewed emphasis 
on the known field-particle duality and long known con- 
cepts as intermittency, dissipative anomaly, direct and 
inverse cascades, usually defined in terms of field charac- 
teristics (the Eulerian description), can now be elegantly 
explained in terms of the statistical properties of particle 
trajectories (the Lagrangian description). In this Letter 
we address the problem of relating the Eulerian proper- 
ties to the Lagrangian ones for active scalar transport. 

The active field a{x^t) advected by an incompressible 
flow v{x^t) is governed by the transport equation 

dta + V • Va = hzA a-\- fa . (1) 

The passive field c(cc, t) evolves in the same flow as 

dtc^v-Vc = f^Ac^fc. (2) 

The two scalars have the same diffusion coefficient 
while fa and fc are independent realizations of the same 
random forcing, with characteristic lengthscale //. The 
difference between active and passive scalars resides in 
their relationship with the flow v: the active field enters 
the velocity dynamics whereas the passive one does not. 
Here, the evolution equation for v is 



dtv + V • = — Vp — AaVa + vAv 



(3) 



Eqs. (1) and (3) are the two-dimensional magnetohydro- 
dynamics equations (see, e.g., Ref. [2]), with the following 
glossary: a is the scalar magnetic potential; the mag- 
netic field is 6 = {—d2a,dia); the coupling — AaVa is 
the Lorentz force (V x 6) x 6; the condition V • v = is 



ensured by the pressure term Vp and u is the kinematic 
viscosity. 

Besides its physical relevance (see e.g. [2]), this ex- 
ample of active scalar is particularly interesting because 
of the conspicuous differences to its passive counterpart. 
Namely, whilst a undergoes an inverse cascade process, 
i.e. it forms structures at increasingly larger scales [3], c 
cascades downscale (see Fig. 1). As a consequence, the 
dissipation of active scalar fluctuations by molecular dif- 
fusivity Ca = /^|Vap vanishes in the limit ^ : no 
dissipative anomaly for the field a. Since the scalar vari- 
ance is injected at a constant rate Fq by the source /a, 
we have that ea{t) = \ J a{x^t)'^dx increases linearly in 
time as ^F^t. On the contrary, passive scalar dissipa- 
tion Cc = Ai:|Vcp equals the input ^Fq and holds c in a 
statistically stationary state (see Fig. 1). 

The probability distribution function (pdf) of a is 
Gaussian, with zero mean and variance Fot (see Fig. 2). 
This is a straightforward consequence of the vanish- 
ing of active scalar dissipation. Indeed, by averag- 
ing over the forcing statistics one obtains: dt{a'^^) = 
n(2n — l)Fo(a^"^~^) (odd moments vanish by symme- 
try) whose solutions are the Gaussian moments: {o?'^) = 
(2n — l)!!(Fot)^. This may be contrasted with the single- 
point pdf of c which is stationary and supergaussian (see 
Fig. 2), as it generically happens for a passive field sus- 
tained by a Gaussian forcing in a rough flow (see, e.g., 
Ref. [1]). 

We now focus on the main purpose of this Letter, that 
is to investigate the statistical properties of active and 
passive scalar in the Lagrangian setting. The fundamen- 
tal property which we will exploit is that the equations 
(1) and (2) can be formally solved in terms of particle 
propagators. Let X{s) be the trajectory of a fluid parti- 
cle transported by the flow v and subject to a molecular 
diffusivity a^, landing at point x at time t. The particle 
moves according to the stochastic differential equation 
dX{s) = v{X{s), s)ds + \/2nd(5{s)^ where (5{s) is a two- 
dimensional brownian motion. The probability density 
p(i/, s\x^t) of finding a particle at point y and time s <t 
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obeys the Kolmogorov equations (see, e.g., [9]): 

-dsp{y, s\x, t)-Vy[v{y, s)p{y, s\x, t)] = f^Ayp{y, s\x, t) , | 

(4) 

dtp{y,s\x,t) + ■ [v{x,t)p{y,s\x,t)]=KA:j,p{y,s\x,t) . 

(5) 

The unusual minus signs in the l.h.s. of (4) are due to the 
fact that particles move backward in time. The solution 
of (1) and (2) can be written in terms of the propagator: 

a{x,t) = ds J dy fa{y,s)p{y,s\x,t) , (6) 

c{x,t) = J dy fc{y,s)p{y,s\x,t) , (7) 

as it can be directly checked by inserting (6) and (7) 
in (1) and (2), respectively, and utilizing (5). It is 
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FIG. 1: Power spectra of active and passive scalar variances 
Ea{k) = 7vk\d{k,t)\'^ and Ec{k) = 7Tk\c{k,t)\'^ . The fluctu- 
ations of a and c injected at the forcing wavelength kf flow 
towards smaller and larger wavenumbers, respectively. At 
k < kf we observe power- law behaviors £;a(A:) - k-'^-o±o-^ ^nd 
Ec{k) - \ while a^t k > kf we find Ea{k) - A:-3-6±o.i 

and Ec(k) ^ A;~^"^^°"^ (for a discussion about the spectra of 
two- and three-dimensional magnetohydrodynamics see, e.g., 
Refs. [3, 4, 5, 6, 7, 8]). In the lower left corner, the fluxes of 
scalar variance Ila,c out of wavenumber k. Negative values in- 
dicate an inverse cascade. In the upper right corner, the total 
scalar variance Ca.c 

(t) = J Ea,c{k,t)dk. The active variance 
ea{t) grows linearly in time whereas ec(t) fluctuates around a 
finite value (see text). The rate of active to passive scalar dis- 
sipation is Ca/cc — 0.005. The data result from the numerical 
integration of eqs. (1-3) by a dealiased pseudo-spectral paral- 
lel code, on a doubly periodic box of size 27r and resolution 
4096^. The forcing terms fa and fc are homogeneous inde- 
pendent Gaussian processes with zero mean and correlation 
(/^(fc, t)/,■(fc^ tO) = [Fo/{2nkf)]6^jS{k + k')6{k - kf)6{t - t') 
where i^j — a,c. The coefficients k and u are chosen to obtain 
a dissipative lengthscale of the order of the smallest resolved 
scales. All fields are set to zero at t = 0, and time is defined 
in units of eddy-turnover time r = If /vrms where // = 2n /kf. 



important to notice that both the active and the pas- 
sive scalar are expressed in terms of the same propa- 
gator p, since a and c are advected by the same ve- 
locity field. Eq. (6) can be written more compactly 
as a{x^t) = {Jq fa{X{s)^ s)ds)x ^ where (. . .)x denotes 
the average over particle trajectories, and similarly for 
eq. (7). 

The interpretation of eqs. (6) and (7) is that the value 
of the scalar field at (x^t) is given by the superposition 
of the input along all trajectories eventually convening at 
that point and time. It is worth remarking that, even in 
the limit of vanishing diffusivity, the particle propagator 
does not collapse onto a single trajectory, as expected for 
incompressible rough flows [1]. The differences between 
active and passive flelds are due to the correlations be- 
tween the propagator p and the forcing which are 
absent, by deflnition of passive transport, between p and 
fc. Indeed, via the coupling term in eq. (3), fa affects 
the velocity fleld and consequently the evolution of the 
propagator p. To further clarify the relationship between 
Lagrangian trajectories and active fleld forcing, it is nec- 
essary to investigate the evolution of the particle propa- 
gator. The problem is that p evolves backward in time, 
according to (4) and the condition p{y^ t\x^ t) = 5{y — x) 
is set at the flnal time t. Contrary to an usual forward 
integration, where v and p can be advanced in paral- 
lel, a brute force backward integration would require the 
storage of velocity conflgurations for the whole lapse of 
integration which is quite unfeasible. To overcome this 
problem we devised a fast and low-memory demanding 
algorithm to integrate numerically eqs. (1), (3) and (4) 
[10]. 

Let us sketch the basic idea of the algorithm. We pro- 
ceed inductively and flrst consider the case t = ti = At, 
i.e. a single step. The flelds a(0) and v{0) and the 
flnal condition p{ti) are given as input (we omit spa- 
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FIG. 2: Pdf's of active (solid) and passive (dashed) scalar 
fields normalized by their standard deviation. The active 
scalar pdf is indistinguishable from a Gaussian (dotted). 
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tial dependences). The fields a(ti) and v{ti) are ob- 
tained by a forward integration wiiose numerical cost is 
denoted Cf. The propagator is then evolved backward, 
at a cost denoted C5, to obtain p{0). The total cost of 
the procedure is Ci = Cf -\- Cb. Consider now the case 
t = In = (2^ — l)At. Our aim is to show that the 
algorithm for 2^ — 1 steps can be split into two succes- 
sive applications with 2^~^ — 1 steps and thus lends to 
a recursive procedure. Integrating forward a(0), v(0) for 
2^"-^ steps give a(tAr-i + At) and i;(tAr-i + At) at a cost 
2^~^Cf. Those fields, together with p{tN)^ are the input 
for the algorithm in the interval \tN-i + At,tAr], which 
gives p{tN-i + At) at a cost Cn-i- A backward step 
is finally needed to obtain p(tAr-i). With this field and 
a(0), v{0) we can now apply the algorithm in the first 
half of the interval [0,tAr-i] to get p{0). The computa- 
tional cost is C5 + Cn-1' Summing up all the costs we 
have the recursive relation Cn = 2^~^c/ + 2Cn-i + C5, 
with solution Cn = y2^c/ + (2^ — 1)05. Comparing 
Cn with the cost of a forward integration of eqs. (1),(3) 
and (5), that is (c/ + C5)(2^ — 1), we have an increase 
by a factor ^ 2(c/+cb) ^' logarithmic in the number 
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FIG. 3: Time runs from bottom to top. First column: time 
evolution of the active scalar field resulting from the numerical 
integration of eqs. (1) and (3). Second column: backward evo- 
lution of the particle propagator according to eq. (4). Third 
column: time evolution of the passive scalar field in the same 
flow. 



of timesteps. As for memory requirements, only the 
couples of fields a{tN — t^), v{tN — t/c), k = 1^ . . . ^ N need 
to be stored - again a logarithmic factor. 

A typical evolution of the propagator is shown in the 
central column of Fig. 3. We can now reconstruct the 
time sequence of the forcing contributions (pa,c{s) = 
J dy fa,c{y^ s)p{y^ s\x^t) which, integrated over 5, gives 
the amplitude of the scalar fields according to (6) and (7). 
As shown in Fig. 4, the time series of (t)a{s) and (l)c{s) are 
markedly different. In the active scalar case, the sequence 
is strongly skewed towards positive values at all times. 
This signals that the trajectories preferentially select re- 
gions where fa has a positive sign, summing up forcing 
contributions to generate a typical variance of a of the or- 
der Fot. Conversely, fc can be positive or negative with 
equal probability on distant trajectories, and the time 
integral in eq. (7) averages out to zero for \s — t\ > r 
(r = If /vrms is the eddy-turnover time). Therefore, typ- 
ically ~ Fqt. The cumulative effect of the correla- 
tion between forcing and propagator is even more evident 
from the relation (j)a{s') ds' = J dy a{y^ s)p{y^ s\x^t) 
derived from (1) and (4). As shown in Fig. 4, the growth 
of the scalar variance is thus related to a strong spatial 
correlation between the propagator and the active field. 
That is further evidenced by comparing the first and the 
second column of Fig. 3: the distribution of particles 
closely follows the distribution of like-sign active scalar. 
This amounts to say that large-scale scalar structures are 
built out of smaller ones which "coalesce" together [3]. 
This situation has to be contrasted with the absence of 
large-scale correlations between the propagator and the 
passive scalar field (second and third column of Fig. 3). 

As remarked previously, the absence of dissipative 




FIG. 4: Top: 0a,c(s) = j dy faAv. s)p{y,s\x,t). The two 
graphs have the same scale on the vertical axis. Here, t = 
32. Bottom: time integrals (j)a{s')ds' (upper curve) and 
(j)c{s') ds' (lower curve). Note the different scale on the 
vertical axis. Recall that (j^a^s') ds' — a{x,t) (and similarly 
for c). 
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anomaly is closely related to the onset of an inverse cas- 
cade. Let us now discuss this issue from a Lagrangian 
viewpoint and consider the squared active field o? . On 
one hand, it can be written as the square of (6). On 
the other hand, one can multiply (1) by 2a to obtain 
the equation dta'^ -\- v • Vo? = + 2a/a — 2ea and 

solve it for eo = in terms of particle trajectories. The 
comparison of the two previous expressions yields 

dsj^ ds' jj fa{y, s)fa{y\ s')p{y, s\x, t)p{y\ s'\x, t) = 

^ dsj^ ds'jj fa{y^ s)fa{y'^ s')p{y', s'; y, s\x, t) , (8) 

where p{y' ^ s'] s|a?, t) = p{y^ s\x^ t)p{y\ s'\y, s) denotes 
the probability that a trajectory ending in (cc,t) were 
in {y,s) and {y'^s'). Integration over y and y' is im- 
plied. To obtain the r.h.s. of (8) consider the expression 
o?{x,t) = 2 J^ds J dy fa{y,s)a{y,s), insert a{y,s) = 
ds' j dy' fa{y' ^s')^ i.e. eq. (6) evaluated at time s, 
and exploit the symmetry under exchange of s and s' to 
switch from a time-ordered form to a time-symmetric one 
to get rid of the factor two. 

Notice as a side remark that by integrating the 
r.h.s. of (8) over cc, exploiting the normalization 
/ dxp{y^s\x^t) = 1 (as follows from(5)), and averag- 
ing over the forcing statistics, one obtains (a^) = F(0)t. 
The result can be generalized to show that {o?'^) = 
(2n — l)!!(F(0)t)^, in agreement with the Eulerian ar- 
gument for the Gaussianity of the single-point pdf of a 
given above. 

Eq. (8) can be recast as {j^ fa{X{s)^ s)ds)\ = 
{[j^ fa{X{s)^ s)ds\^)x^ leading to the conclusion that 

fa{X{s), s)ds is a non-random variable over the en- 
semble of trajectories. This result has a simple interpre- 
tation: the absence of dissipative anomaly (i.e. Ca = 0) 
is equivalent to the property that along any of the in- 
finite trajectories X{s) ending in (x,t) the quantity 
Jq /a(X(s), s)ds is exactly the same, and equals a{x, t). 
Therefore, a single trajectory suffices to obtain the value 
of a{x^ t) contrary to the passive case, where different tra- 
jectories contribute disparate values of fc{X{s)^ s)ds^ 
with a typical spread €ct. In this case, only the average 
over all trajectories yields the correct value of c{x^ t). 

The previous property extends to multiple points. Pro- 
ceeding as for (8), one can show that m trajectories are 
enough to calculate the product of arbitrary powers of a 
at m different points. In particular, structure functions 
((a(r) — a(0))^) for any order n can be shown to involve 
only two trajectories. This should be contrasted with the 
passive scalar, where the number of trajectories increases 
with n and this is at the core of anomalous scaling of the 
field [1]. Indeed, we do not observe any anomalous scal- 
ing for the structure functions of the active field a and 



the pdf's of its increments are self-similar. Note that 
the same phenomenon and a lack of dissipative anomaly 
were also found for a passive scalar in a compressible 
flow [11]. In that case, the ensemble of trajectories col- 
lapses onto a unique path, fulfllling in the simplest way 
the constraint (8) and its multi-point analogs. Here, even 
though the trajectories do not collapse, the constraints 
are satisfied due to the subtle correlation between forcing 
and trajectories peculiar to the active case. 

We conclude with the following open issue. It is known 
that passive correlation functions are controlled by the 
asymptotically dominant terms in particle propagators 
(see, e.g., Ref. [1]). Is this also true for active scalars ? 
The correlations between propagator and input in (6) 
suggest this be not the general rule. For some specific 
systems, those correlations might however be such that 
the asymptotics of the propagator still stands out at large 
times [12, 13]. Clarifying the dynamical conditions con- 
trolling this phenomenon will be the subject of future 
work. 
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